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Abstract 



> 

CN ■ We present general algorithms to convert scattering data of linear and area detectors 

recorded in various scattering geometries to reciprocal space coordinates. The pre- 
sented algorithms work for any goniometer configuration including popular four-circle, 
six-circle and kappa goniometers. We avoid the use of commonly employed approxi- 
mations and therefore provide algorithms which work also for large detectors at small 



sample detector distances. A recipe for determining the necessary detector parameters 
C3 ■ 

including mostly ignored misalignments is given. The algorithms are implemented in 
a freely available open-source package. 



1. Introduction 

Elastic X-ray scattering is a very widely applied technique to study the structure of 
materials ranging from single crystals, powders and other forms of hard condensed 

PREPRINT: Journal of Applied Crystallography A Journal of the International Union of Crystallography 



2 
matter to biologic tissues and organic molecules. Crystalline as well as non-crystalline 
matter like liquids and amorphous materials can be studied by techniques like X-ray 
diffraction and reflectometry. A variety of approaches exist to analyze the scattering 
data. Some quantities can be directly determined from the measurements (e. g. lat- 
tice parameters (Bond, 1960; Fewster, 1999), layer thicknesses (Warren, 1969; Pietsch 
et al., 2004)). Other types of analysis involve simulation of the scattering signal 
to determine strain and material composition (Stangl et al., 2004; Wintersberger 
et al., 2010), or the model- free determination of real space structure using phase 
retrieval algorithms (Fienup, 1982; Miao et al., 1998; Diaz et al., 2009; Minkevich 
et al., 2011). All of those approaches have in common that the analysis is most of 
the time performed in reciprocal space, and hence require a conversion of experimen- 
tally measured data into reciprocal space. While the particular analysis steps differ 
for different experiments, the conversion itself is a common step, which needs to be 
performed for a lot of different techniques. This has been treated very well for the case 
of point detectors (Busing & Levy, 1967; Lohmeier & Vlieg, 1993; You, 1999; Bunk & 
Nielsen, 2004). For ID and 2D detectors, which are used more and more frequently, 
several issues related to the detector geometry and detector alignment complicate a 
correct conversion. 

For some application like protein crystallography or powder diffraction, experi- 
mental as well as analysis schemes are standardized. Examples are software pack- 
ages to extract peak positions and intensities in protein-crystallography(Leslie, 2006) 
and powder diffraction (Lutterotti et al., 1999; Rodriguez- Carvajal, 2001; Hammers- 
ley, 2013), or the pdb protein data bank format and cif crystallographic information 
file format established by the International Union of Crystallography for exchange of 
structure files and experimental data. 

For most other cases, only very specialized solutions to particular experiments exist, 
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each containing solutions for some aspects important for the respective case, e. g. the 
experimental control software spec (spec, 2013) used at various synchrotron sources 
is able to do reciprocal space conversion for several goniometer geometries, however 
only for point detectors. Commercial software supplied with several diffractometers is 
optimized for certain geometry /detector combinations. 

A particular problem of ID and 2D detectors are misalignments with respect to 
the ideal case: At zero detector angle, the line or plane of the detector should ideally 
be perpendicular to the incident X-ray beam. In practice, this will not be the case, 
even if deviations will usually be very small. For 2D-detectors, in addition a rotation 
of the detector around the incident beam direction can occur. In most cases, these 
misalignments will be unintentional and rather small, they are difficult to measure 
and hence often neglected. However, considering these effects is actually important to 
obtain correct and accurate reciprocal space coordinates. 

We present a general applicable algorithm for the conversion of experimental data 
recorded with point, linear and area detectors, and for any diffractometer with arbi- 
trary number and sequence of axes. For this purpose, we generalize the algorithms 
presented in several papers (Busing & Levy, 1967; Lohmeier & Vlieg, 1993; You, 1999; 
Bunk & Nielsen, 2004) to arbitrary goniometer geometries without approximations: 
we consider the fact that most ID and 2D detectors are flat and hence the relation 
between channel number and scattering angle is non-linear. Furthermore, we provide 
recipes how to determine the necessary detector parameters from a set of simple scans 
around the primary beam. These scans also enable the user to determine the men- 
tioned and several other misalignment parameters. 

To keep our solution as general as possible, we have implemented it within a freely 
available software package xrayutilities (Kriegner & Wintersberger, 2013). This gen- 
eralized algorithm is also particularly useful for automatized tool-chains as planned 
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4 
by several synchrotron facilities right now (Passerelle, 2013). In addition to the recip- 
rocal space conversion described in this article, xrayutilities includes routines to read 
various data formats and methods to calculate experimental parameters from mate- 
rial properties. More information on those parts can be found on the xrayutilities 
web page (Kriegner & Wintersberger, 2013). This article focuses on the reciprocal 
space conversion part: After the introduction of the applied algorithms following the 
approach of You, 1999 we show the extension for linear and area detectors and explain 
the use of our algorithms and how necessary parameters can be determined. We demon- 
strate the application of our approach for few selected examples of both laboratory 
diffractometers and synchrotron beamlines. In the appendix we discuss one complete 
example, including the particular entries into the xrayutilities package to define the 
diffractometer geometry, correctly initialize a 2D detector setup, and convert a 2D 
detector frame into reciprocal space. 

2. Angular to reciprocal space conversion 

Conversion of angular coordinates to reciprocal space can be tedious since one needs 
special equations for every diffractometer/detector geometry. For several diffraction 
geometries explicit formulas are given for example in Pietsch, 2004. However the con- 
version can be performed in a general way as long as the information about the 
goniometer geometry is available together with the experimental angles. The algorithm 
presented below therefore needs not just the experimental angles as input parameters 
but also the description of the goniometer. To work for arbitrary goniometers the 
physical order of the cradles, i. e. how they are mounted on each other, as well as the 
orientation of the rotation axes are needed. An example is given below. To unambigu- 
ously specify the rotation axes of the goniometer circles, a reference coordinate system 
is fixed to the laboratory frame. It is useful to choose this coordinate system in a way 
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5 
that the primary X-ray beam propagates along one of the coordinate axes. 

The reciprocal space coordinates we want to know for each measured point are 
coordinates of the momentum transfer q = kf — fcj , where ki is the wave- vector of the 
incident X-ray beam, kf the wave-vector of the scattered beam towards a particular 
detector (pixel) position. The coordinates of q we are interested in are those in a 
coordinate system fixed to the sample under investigation. The coordinates of fcj in 
the laboratory system are fixed by our choice of the coordinate system. The coordinates 
of kf are given by the angle positions of the detector arm. To describe q in the sample 
coordinate system, we also need all goniometer angles changing the sample orientation. 

A minimal two dimensional example illustrating the definition of our laboratory 
coordinate system (blue) is shown in Fig. 1. The coordinate system attached to the 
sample is indicated in green in Fig. 1. So what we have to deal with are coordinate 
transformations between the different involved coordinate systems. Those coordinate 
transformations can be written as matrix equations. We will reproduce some of the 
essential equations of earlier papers (Busing & Levy, 1967; You, 1999) so that reader 
can follow the further generalization for one and two dimensional detectors. We con- 
sider just a point detector for the moment, the generalization for one or two dimen- 
sional detectors is shown in Sec. 3 below. 

The sample coordinate system we have been talking about is actually the coordi- 
nate system attached to the innermost goniometer circle (because this is what the 
goniometer angles describe). Of course, the sample can in addition have a certain ori- 
entation with respect to this circle. This is most evident in the case of a crystalline 
sample, where the directions of the reciprocal space of the crystal have a certain orien- 
tation with respect to the sample holder. This and the coordinates of the momentum 
transfer in the reciprocal space of the crystal are described below by matrices U_ and 
22, respectively. In the following, we describe this "most complicated" case of a sin- 
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6 
gle crystalline sample. For amorphous or powder samples etc., U_ and/or B can be 
replaced by identity matrices. 

In the crystal lattice a momentum transfer might be described by the column vector 

3 

h c = Y J hjbj = Bh (1) 

3=1 

with bj the reciprocal lattice vectors or the matrix B formed from those vectors and hj 
generally referred to as h, k, I. In case h, k, I are integer values they are called Miller 
indices. An explicit representation of £ is for example given in (Busing &; Levy, 1967). 
To connect the momentum transfer in the crystal with the momentum transfer in 
the coordinate system attached to the sample goniometer the orientation matrix £/ is 
introduced. 

h s = Uh c = RBh (2) 

This momentum transfer converted to the laboratory frame is equal to q. The con- 
version of the momentum transfer is achieved by another coordinate transformation 
(described by matrix S), which solely depends on the sample orientation and therefore 
the goniometer angles which move the sample. 

q = k f -ki = £h s (3) 

In the case of integer h, k, I, the former equation basically is the Laue equation. 

Also any detector rotation can be expressed as coordinate transformation described 
by i2 and therefore the exit wave vector is given by 

k f = Rh (4) 

Combining Eqs. 2 to 4 one obtains the momentum transfer in the crystal coordinate 
system 

hc=(§ = Iiy 1 (R-l)k i (5) 
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7 
with the identity matrix 1. Note that all of the used coordinate transformations 
are transformations between two orthogonal coordinate systems, and therefore yield 
orthogonal matrices, which are invertible. The matrix 22, which is formed from the 
reciprocal space vectors bj, which are linear independent is also invertible. Since the 
above conversion involves a matrix inversion this is important for the algorithm to 
work. 

The rotation matrices £ and 22 can be deduced from the description of the goniome- 
ter by multiplying the rotation matrices from each circle starting with the outer most 
circle. 

22, S. = (outer most) (inner most) (6) 

For that purpose the goniometer rotation axis needs to be defined in the laboratory 
coordinate system for the case when all circles are set to zero. It is therefore useful to 
choose the coordinate system in a way, which allows this description to be as simple as 
possible. Keep in mind that later also the detector directions need to be determined 
in this coordinate system. In addition to the direction of the rotation axis also the 
rotation sense needs to be described. We use the mathematical definition of rotation 
sense. For most goniometers (except for special geometries like K-goniometer, treated 
separately below) the rotation axis point along primitive directions. When looking on 
the rotation from the positive side of the corresponding direction, clockwise rotation 
is left-handed, i.e., negative, and an anticlockwise rotation is right-handed (positive). 
For example, we will call a clockwise/left-handed/negative rotation of angle a around 
the x-axis a 'x-'-rotation, described by the following rotation matrix: 



1 











cos a 


sin a 





— sin a 


cos a 



R (a) = cos a sin a (7) 
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2.1. example of a goniometer definition 

To elucidate the definition of a goniometer we use the goniometer shown in Fig. 2. 
The goniometer has a 3S+2D geometry, which means it offers three degrees of rota- 
tion for the sample and two independent degrees of freedom for the detector. The 
goniometer axes are specified by their axis direction and rotation sense. The coordi- 
nate system is chosen to have x pointing downstream along the primary beam, z is 
pointing upwards and y pointing backwards to have a right-handed reference frame. In 
xrayutilities we describe each rotation axis with one character giving the axis direction 
(either x, y, or z) and the rotation sense (either + or -). This description needs to 
be supplied for the sample and detector circles for the case where all axis are at zero 
position starting with the outermost circle. The goniometer in Fig. 2 has three sample 
circles (fx, x, c / 3 ) with the indicated rotation directions. The outermost angle fi turns 
clockwise around the z-coordinate axis and thus is described by 'z-'. The complete 
sample goniometer is described by the following set of rotations: ('z-','x-','y+'). For 
the detector circles turning around the z and y-directions, the corresponding defini- 
tion is ('z-','y-'). A f un example of how to insert these definitions into xrayutilities is 
given in the appendix. 

2.2. special rotation directions (Kappa goniometer) 

In addition to rotations around axis of the coordinate system an often used goniome- 
ter geometry is the so called K-geometry (Poot, 1972; Thorkildsen et al., 1999), in 
which one of the rotation axis has an angle of typically 45 to 60° with one of the other 
axis. In xrayutilities we define such an axis by k+ or k-. The plane of the /{-rotation 
axis as well as its angle with respect to a reference direction are specified in a configu- 
ration file by the options kappa_plane and kappa_angle. The rotation matrix needed 
in matrix S. can be determined easily using an general equation given in the appendix. 
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3. Angular to momentum space conversion for ID and 2D detectors 

Equation 5 describes the conversion from angular coordinates of a general goniometer 
to reciprocal space for a point detector only. The generalization for a linear or area 
detector requires information about the pixel size and distance, and the direction 
into which a pixel row extends. Usually in the data files only one angular position is 
stored for every data point recorded with a multidimensional detector. This angular 
coordinate corresponds to the position of the so called center channel/pixel (no), which 
is the pixel hit by the primary beam when all the angles are set to zero. For all other 
detector pixels their position needs to be determined. This is easiest for the case of a 
curved one dimensional detector in which every detector channel or pixel covers the 
same detector angle segment. Every detector pixel (identified by the channel number 
n) corresponds to a detector angle 26 given by 

26>(n) = 26> + (n - n )/N (8) 

where TV is the number of channels per unit of angle of the detector circle. When the 
detector angle is expressed in degree N equals the number of channels per degree of 
rotation. 

Most modern detectors are, however, straight as shown in Fig. 3 and not curved and 
therefore Eq. 8 is not generally applicable. Only in the limit of a large sample detector 
distance the curved and straight detectors become indistinguishable. Nevertheless the 
channel per degree approximation is frequently used in practise. In xrayutilities one- 
dimensional detectors are treated as straight detectors and Eq. 5 is adjusted accord- 
ingly. For each detector pixel n, the corresponding direction of a scattered beam hitting 
this pixel is calculated, replacing Eq. 4 by 



kf(n) = \kj\ D (kj + d(n — no)w/L) /\\ki + d(n — no)w/L\\, (9) 

where w and L are the size of a detector pixel and distance from sample to detector as 
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shown in Fig. 3. The direction in which the detector channel number increases is given 
by d. A 'hat' on a vector indicates a unit vector. The fraction w/L is in the case of 
large sample to detector distance equal to 1/N, where N is is the number of channels 
per radians, and Eq. 9 effectively simplifies to the form of Eq. 8. For a two dimensional 
detector with channel directions d\ and d,2 we can write an equivalent equation for 
the exiting wave vector of channel (ni,ri2) including the width of the detector pixels 
wi t 2 as 

fe/(ni,n 2 ) = \ki\R[ki + d\{n\ - n\ ')w\/L + d 2 {ri2 - n 2 >)w2/L) /|| ••■ ||, (10) 

Using the description of the detector in real space we do therefore avoid the " channel 
per degree" approximation, which implicitly assumes a detector is curved and therefore 
does not work for small sample detector distances with large detectors. Inserting Eqs. 9 
and 10 into Eq. 2 and 3 yields general equations for the reciprocal space conversion of 
linear and area detectors. The possible misalignments are included in those equations 
via the pixel directions d for the ID and d\p for the 2D detectors. 

3.1. Detector parameters of ID detectors 

For the conversion algorithms described above several detector parameters are 
needed. Among them the detector distance L and width of one detector channel w. 
Since neither the detector distance nor the width of one channel can be easily mea- 
sured very accurately, we use the fact that only their ratio is needed, which can be 
determined from an angular scan through the primary beam. Assume a linear detec- 
tor is mounted along the direction in which the detector arm of the used instrument 
moves. Scanning the detector angle will move the primary beam over the detector, 
by modeling this movement we are able to determine the needed quantities. There- 
fore assume a linear detector mounted at a distance L like the one shown in Fig. 3. 
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Neglecting for the moment a possible detector tilt, i. e. when a detector is not mounted 
perfectly perpendicular to the X-ray beam, the detector channel at which the detector 
is hit for a detector arm angle 29 is given by 

n (20) = — tan26» + n . (11) 

w 

If a detector tilt ft as shown in Fig. 3 is included above equation needs to be modified 

and one finds 

, ~ ^ L sin 29 

n(2$,P) = — — - + no. 12 

it; cos (29 — p) 

By fitting one of the two models one can find the necessary detector parameters needed 

for the reciprocal space conversion of a linear detector. For this purpose a scan through 

the primary beam should be performed with the linear detector and the detector 

spectra should be saved at every position. Since not only the slope is determined 

(which only needs two spectra), it is required to acquire several spectra, we suggest 

typically > 10. For the determination of the parameters two functions are provided 

in xrayutilities. One of them automatically processes the spectra of a linear detector 

acquired during a scan through the primary beam and determines the position of the 

primary beam in every spectrum by fitting a Gaussian peak. From this fitting the 

position of the primary beam is found with sub-pixel precision. The second function 

needs the user to determine the channel number of the primary beam and supply it to 

the function, which should be used only in cases the first function is not applicable. 

Calling one of the two functions will produce a plot as the one shown in Fig. 4, which 

shows the channel number at which the beam was hitting the detector together with 

the variation expected from the model(s). A second subplot shows the comparison of 

model and measured data with the linear trend subtracted, for two different detector 

distances of 380 mm and 250 mm for a straight linear detector with a pixel size of 

50 //m. The different distance shows up as different slope in the upper subplot. When 
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the linear trend is subtracted the non-linearity due to the trigonometric functions in 
Eqs. 11 and 12 can be seen. The fit is also sensitive to a tilt of the detector from the 
direction perpendicular to the primary beam as can be seen from the comparison of 
the model without tilt (Eq. 11, dashed line) and model with tilt (Eq. 12, full line). 
For the measurements shown in Fig. 4 an artificial tilt of 0.3 deg was introduced, 
which was also determined by the fit. Note that a tilt of the detector is not only 
a result of a not perfectly mounted detector, but can also result from the fact that 
experimentalists choose to use a center channel, which does not correspond to the true 
center of the detector and do this by redefining the zero position of the detector arm 
rotation, which introduces a tilting of the detector by exactly the angle by which the 
detector arm angle changed! Using a translation of the detector perpendicular to the 
beam (mounted on top of the detector arm) would prevent such a tilt, is, however, 
very often not available. 

Basically this shows the necessity of non-linear models instead of the simpler lin- 
ear fitting, which is a reasonable approximation only in the case of a large sample 
to detector distance. If a linear or area detector covers an angular range of > 8° it 
gets absolutely necessary to use the above treatment since deviations would exceed one 
channel for typical channel sizes of around 50 /im. Using large linear detectors at mod- 
erate sample detector distances this limit is frequently reached especially in modern 
laboratory diffractometers. To further highlight the necessity of the non-linear models 
we show an example of a x-ray diffraction reciprocal space map measurement of the 
Si (331) Bragg peak of a Si(lll) oriented single crystal in Fig. 5. The measurement 
was performed with the above mentioned linear detector at a distance of 250 mm at a 
laboratory diffractometer with CuKal radiation produced by a Ge(220) channel cut 
monochromator. The same measurement was repeated three times, however, different 
parts of the detector were used for the detection of the diffracted signal. We compare 
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the reciprocal space maps obtained using our exact conversion formalism in compari- 
son with the "channel per degree" approximation which assumes a curved detector like 
described by Eq. 8. Using this approximation (panels (d) to (f) in Fig. 5) we find that 
only the measurement using the central part of the detector gives the correct peak 
position in reciprocal space/lattice parameter. The measurements performed with the 
detector offset to higher or lower angles would result in a lattice parameter wrong by 
approximately 0.04%, which is far beyond the usual sensitivity of the experimental 
setup, which is said to be < 1 x 10~ 4 and can be increased further as for example 
outlined in (Fewster, 1999). Using the accurate conversion no such shifts are observed 
and the three measurements (panels (a) to (c)) are indeed equivalent. 

3.2. Detector parameters of 2D detectors 

For 2D detectors similar determination of the wx^jD ratio can be performed if the 
detector rotation and other misalignments (see below) are negligible, by decomposing 
the problem into two separate one dimensional problems. In case the detector is rotated 
around the primary beam the problem can no longer be decomposed. In particular one 
more problem specific to two dimensional detectors arises. Very often the true zero 
of the outer detector arm rotation is not known. For the inner detector rotation this 
does not imply any particular problem since it shows up only as additional tilt, which 
can be easily corrected. However, an offset in the outer detector rotation implies a 
rotation of the rotation axis of the inner detector rotation. In case the outer motor 
is not at the true zero the inner rotation will no longer rotate perpendicular to the 
primary beam as indicated by the two rotation planes shown in Fig. 2. In case the 
outer motor would be at its true zero the detector would rotate in the blue plane, due 
to an offset in the outer rotation the detector instead rotates in the red plane. This 
means the offset of the outer motor needs to be determined from alignment scans as 
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well. 

In general one needs to determine 8 parameters for a 2D detector: the center chan- 
nels (n\ , n 2 ), the ratios w\/L and W2/L, and directions of the vectors d\p, which 
are specified by the tilt, tilt direction and rotation of the detector around the pri- 
mary beam and the outer angle offset described above. The untilted directions of the 
vectors d\ t 2 can be determined by knowledge of the primary beam and detector rota- 
tion directions unambiguously and do therefore not need to be included in the fit. 
Similar as for the one dimensional detector those parameters can be determined from 
scans through the primary beam. It is necessary to use at least two scans, one with 
the inner and one with outer detector arm rotation with sufficient points around the 
primary beam and acquire a detector image at several positions. From the primary 
beam position observed in those images the detector directions and other necessary 
parameters can be determined by a fitting routine. As quality criterion of the fit the 
average ||qr|| -value of the pixel at which the primary beam is observed is used. This 
position is calculated using Eqs. 3 and 10. Since it is the primary beam we observe in 
all the images, this ||q||-value should be zero when the correct detector parameters are 
found, xrayutilities provides a function which takes the detector images and performs 
a fit of the 8 described parameters. Since several of those parameters, e. g. the offset of 
the outer motor with one of the center channels, are coupled with each other, the fit 
is performed in a way that it starts from several starting parameters to find the global 
minimum in the parameter space. An example of such a fit is shown in Fig. 6. The 
detector parameters have been determined from two perpendicular scans through the 
primary beam using a Maxipix-detector with effectively 516 times 516 pixels at beam- 
line ID01/ESRF (Microdiffraction imaging beamline ID01 - ESRF Grenoble, 2013). 
One scan was performed using the outer detector arm motor (moving horizontal) and 
the other one with the inner detector circle (moving vertical). In each scan we used 70 
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points, which means that in total 140 images were used. Note that using considerable 
less number of points does reduce the quality of the fit until no reliable statement 
can anymore be made about the misalignment parameters. It is therefore suggested to 
use a comparable number of points as in this example. Furthermore the images were 
manually selected to use just those with the full primary beam in active area of the 
detector. These images were fed into the algorithm described above and the detector 
rotation, tilt, tilt azimuth and outer angle offset along with the center channels and 
detector pixel size were determined. This determination is shown in Fig. 6, were the 
average ||g|| -deviation of the detector pixels positions associated with the primary 
beam in the performed scans is shown. The value of this deviation is approximately 
the offset of the absolute value of the momentum transfer in reciprocal space. This 
deviation shows a clear minimum in all the 8 parameters. These 8 parameters are: 
1/2. cchl,2 are the center channels of the detector in both directions. This is the 
pixel number where the primary beam is hitting the detector at the true zero position 
of the detector arm (including the outer angle offset). 3/4. pwidthl,2 are the pixel 
widths of the pixels in the two detector directions. The unit of these values in the plot 
is the size of the pixels in jjlvh assuming a detector distance of 1 m. This corresponds 
to the parameter w/L from above. If the pixel size is known the detector distance 
can be calculated or vice versa. 5. tiltazimuth is an angle giving the direction of the 
detector tilt. An value of 90 and 270 degrees corresponds to a tilt rotation around the 
first pixel direction and and 180 around the second pixel direction, respectively. 

6. tilt is the tilt angle of the detector plane around the direction given by the 
tiltazimuth, since the tilt-azimuth runs from to 360 only positive tilts are used. 
7. detrot is the detector rotation around the primary beam direction in degrees. 8. 
outerangle offset is the offset of the outer detector arm rotation in degrees. 

For the determination of those parameters fits with various different starting param- 
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eters are performed. This is absolutely mandatory since several of the parameters are 
correlated, and therefore a single fit would not necessarily find the global optimum. 
One correlation, which can be easily imagined is the correlation of one of the center 
channels with the outer rotation offset. Since the detector orientation is not given to 
the script and automatically determined from the given scan data it is not apriori 
clear if this is center channel 1 or 2. In Fig. 6 this correlation can be seen in the plot 
showing center channel 2 (cch2) and the outer angle offset. The cloud of points from 
center channel 2 is rather broad since this value changes when the outer angle offset 
is changed. In fact the coloring of the cloud of points of center channel 2 reveals that 
it is only mirrored with respect to the one of the outer angle offset. A not so intuitive 
correlation is the one of the detector tilt with the outer angle offset. Imagine a not 
tilted detector, which is however offset in the outer angle. Such an offset will effectively 
tilt the detector by exactly the offset in the outer angle with a tilt azimuth of 90 or 
270 deg. In case the detector is already slightly tilted in an arbitrary direction without 
outer angle offset, the tilt due to an outer angle offset will overlay with the initial tilt 
and influence the tilt azimuth and tilt in a non-trivial manner. The optimal set of 
parameters for the area detector of beamline ID01 was determined from the point 
with the lowest lowest error (global minimum in the parameter space) of 2.70 x 10~ 9 . 
To elucidate the benefit of considering those "misalignment" parameter, which are 
usually omitted, we also give the obtained errors when several of the quantities are 
fixed. When only the center channels and pixel widths are fitted we obtain an error of 
1.68 x 10~ 6 , which is orders of magnitude higher than our optimal error. In the present 
example the most important parameter is the detector rotation which brings the error 
down to 3.65 x 10~ 9 . The second most important parameter is the outer angle offset, 
which brings down the error further to 2.81 x 10~ 9 . If only the tilt (tilt azimuth and 
tilt angle) or the outer angle offset are fitted without the detector rotation, the error 
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can only be reduced by less than 10% from the value without considering any mis- 
alignment. Only when all parameters are considered in the fit the error can be reduced 
to the optimum. This clearly shows that all those parameters should be considered to 
obtain the correct reciprocal space positions of the full area detector. 

It should be noted that some of the parameters like the offset of the outer detector 
arm rotation and detector tilt can only be determined when the detector is spanning 
a certain angular range, the user is therefore urged to check the resulting plots in 
order to see if the parameters are reasonable. In cases where the detector distance is 
large the outer angle offset can only be determined with large error and it might be 
better to fix the corresponding parameter during fitting. Also the in cases when the 
detector tilt is small the tilt azimuth will not have a clear minimum and is therefore a 
indeterminable parameter. The script also outputs the code needed for initialization of 
the area detector, which considers the determined tilts and rotations as shown in the 
appendix B. These detector tilts of linear and area detectors are then included in the 
reciprocal space conversion. When the detector pixel position is calculated as described 
in Eq. 9 and 10 the detector direction d or di ; 2 need to be rotated accordingly before 
the detector position is calculated. 

In case no detector rotations are available, e. g. when using an area detector for 
powder diffraction we refer to the fit2d software (Hammersley, 2013) for determining 
the necessary parameters. 

4. Conclusions 

In conclusion we present algorithms for reciprocal space conversion for X-ray diffrac- 
tion experiments. We generalize the equations given by (Busing &; Levy, 1967; You, 
1999) for the use of linear and area detectors. Using our approach we can convert angles 
from arbitrary goniometers to reciprocal space coordinates. For the conversion of linear 
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and area detectors several detector parameters, including all possible misalignments 
are needed. Recipes were presented how those parameters can be determined from 
alignment scans for both linear and area detectors. The software package including 
the presented algorithms is available at http://xrayutilities.sourceforge.net. 
The algorithms were shown to determine the detector parameters of a ID and 2D 
detectors, which are otherwise not determinable. 



Appendix A 
Implementation details 



The software package is available from http://xrayutilities.sourceforge.net 
and is mainly coded in the popular scripting language Python™ with some perfor- 
mance critical parts written internally in the C programming language. The user needs 
to use only the Python package. Installation is possible on all major platforms (Win- 
dows, Mac, Linux/Unix). In addition to the reciprocal space conversion the package 
offers several functions to read data from spec-files, spectra-Hies and from laboratory 
diffractometers of Panalytical (xrdml) and Seifert (nja) as well as several CCD-data 
files (edf,tiff) and other function, which are described in the documentation found at 
the web-page. 

The reciprocal space conversions performance critical part, which implements Eq. 5 
and equivalents for linear and area detectors is coded in C. As can be seen from Eq. 5 
it mainly consists of 3 x 3 matrix operations. For the conversion of angular positions 
the matrices of the sample rotations need to be set up, for which the rotation matrices 
as the one shown in Eq. 7 are used. These matrices are multiplied in the order as 
given in Eq. 6. Furthermore, the conversion involves one matrix inversion, which is 
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performed using the adjugate matrix formula. For matrix A this means the inverse is 

given by 

adj A 

A' 1 = -3 t • (13) 

= detA 

Depending on the number of sample (n s ) and detector circles {rid) to be consid- 
ered in the conversion for a point detector n s + rid + 2 matrix multiplications, the 
matrix inversion, and one matrix-vector multiplication as well as one matrix-matrix 
subtraction need to be performed for every data point. 

However, for linear and area detectors not for every detector pixel all the conversion 
needs to be done. For one spectrum or image all the sample and detector angles are 
considered to be constant and therefore the n s + rid + 1 matrix multiplications as 
well as the matrix inversion need to be performed only once. Only two matrix-vector 
multiplications and some vector arithmetic to set up the pixel position need to be per- 
formed for every detector channel/pixel which enables a fast conversion also for large 
area detectors. Furthermore, the conversion for multiple data points is independent 
and can therefore be easily performed in parallel on modern multiprocessor computers. 

Appendix B 
example script for the reciprocal space conversion 



In the following an example script for the conversion of an image from an area 
detector recorded using the goniometer shown in Fig. 2 is given. This goniometer 
has the sample angles /x, x an d <f> and detector angles v and 5. The rotation axis 
of this goniometer are as given in the main text (sample circles: ( 'z-' , 'x-' , 'y+O 
and detector circles: ('z-' , 'y+')) and the primary beam propagates along positive 
x-direction. The script is written in the Python™ script language. 

import xrayutilities as xu 
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Nav = (1,1) # number of pixel to average over (reduce amount of data) 
region_of .interest = (100,500,100,500) # region of interest on detector 
en = 9000 # xray energy in eV 

qconv = xu.QConversion(( ' z-' , 'x-' , 'y+' ),( ' z-' , ' y-' ), (1,0,0) ) 

# 3S+2D goniometer: sample rotations: mu,chi,phi: ( 'z-' , 'x-' , 'y+') ; 

# detector rotations: nu, delta: ('z-'.'y-') 

# primary beam direction (1,0,0) 

Si = xu. materials .Si 

hxrd = xu. HXRD (Si. Q( 1,0,0) , Si. Q (0,1,0) ,en=en,qconv=qconv) 

hxrd.Ang2Q.init_area(>z-> , 'y+' ,cchl=300 . 11 ,cch2=320.78,Nchl=516,Nch2=516, 
pwidthl=l . 6639e-04 , pwidth2=l . 6630e-04 , distance=l . , detrot=-0 . 749 , 
tiltazimuth=3 . , tilt=0 . 448 ,Nav=Nav ,roi=region_of .interest) 

outerangle_of f set = -0.643 

# experimental goniometer angles 
mu,chi,phi = (20,0,0) 

nu, delta = (40,0) 

# conversion to Q 

qx,qy ,qz = hxrd. Ang2Q .area(mu, chi,phi,nu,delta,delta=(0,0,0,outerangle_of f set ,0) ) 

# conversion to (h,k,l) 

h,k,l = hxrd. Ang2HKL (mu , chi , phi, nu, delta, delta=(0,0 ,0,outerangle_of f set ,0) , 
mat=Si,dettype='area' ) 



All the detector parameters (like center channels, number of pixels, size of the 
pixels and the distance as well as possible tilts and rotations of the detector) need to 
be given to the init_area-function. This parameters are the result of a fit like the 
one shown in Fig. 6. The conversion is then performed by calling the area-method of 
an experiment class object, where also the outerangle offset needs do be considered. 
An experiment object holds information about the experiment like orientation of the 
crystal on top of innermost sample circle, the goniometer geometry and the X-ray 
energy. Here the object is initialized for an experiment performed using X-rays with 
an energy of 9000 eV and silicon crystal with the [100] direction along the primary 
beam at zero position and [010] along the z-axis. The second argument [010] gives the 
reference direction perpendicular to the primary beam in the plane spanned by the 
innermost detector direction. Dummy values for the reference directions can be used 
when the Ang2HKL function is not used! 
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After execution of this script the variables qx, qy and qz hold the reciprocal space 
position of the detector pixels specified by the region of interest variable. In the case 
above those variables will have a shape of 400x400, holding the (/-position of the pixels 
100 to 500 in both detector directions. The variable Nav could be used to reduce the 
amount of data by block averaging. Using Nav = (2,2) effectively quadruples the area 
covered by one pixel and therefore reduces number of returned (/-positions to 200x200 
per coordinate. 

Furthermore the use of the algorithms to directly convert angles to reciprocal space 
indices h, k and / is shown. A function named Ang2HKL, which is called together with 
the orientation matrix to directly convert the experimental angles to the indices h, k, I 
instead of the momentum transfer in units of A" 1 . Instead of specifying the orientation 
matrix one can also use the crystallographic orientation of the sample to determine the 
matrix £/, which is the default if no ^-matrix is explicitly given. Special algorithms 
also exist for determining the matrix B_ from unit cell vectors. In the example the 
matrix 22 is determined directly from the built in Material properties of silicon. In the 
present case (horizontal diffraction) the momentum transfer for the area detector has 
mainly a component along negative y-direction, whereas the h, k, I components have 
there largest component along I. 

Experimental angles do of course not need to be entered manually. They can be 
read from several different data formats as mentioned above. Details are given in the 
documentation found online. 



Appendix C 
rotation matrix for arbitrary axis direction 



Several rotations need to be performed around none primitive directions, e. g. in the 
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case of K-goniometers or also when calculating detector tilts. In such case the rotation 
axis as well as the angle of rotation are known and rotation matrices are generated 
using an formula given in (Lang h, Pucker, 1998). The components or the rotation 
matrix (J£)ij = r^ for rotation of angle a around the axis e are given by 



Tij — 6 



e j + (<% ~ e i e j) cos a ~ e ijkek sin a (14) 



with the Kronecker delta 5ij and the Levi-Civita symbol e^. 
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detector n 




oratory 



Fig. 1. Two dimensional sketch of the scattering process and the coordinate systems 
attached to the laboratory (xl,Ul) and sample (x s ,y s ). Shown are the incidence and 
exit wave vectors kij as well as the momentum transfer of the scattering process 
q. The laboratory coordinate system is chosen to have the x direction to coincide 
with the primary beam direction. 
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Fig. 2. sketch of a goniometer with three sample axis, and two detector rotations. 
A red and a blue plane indicate the detector rotation planes of the inner detector 
rotation for two different positions of the outer detector arm rotation. 
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Fig. 3. Sketch of a linear detector mounted at distance L from the center of rotation 
of the goniometer. In (a) the detector direction specifying the direction along which 
higher channel numbers are found is given by d. Also indicated is the width of 
one channel w and the center channel no, which is the channel where the primary 
beam is centered at zero detector angle. Panel (b) shows the possible detector tilt 
f3 (misalignment), resulting if a detector is not mounted perfectly perpendicular to 
the X-ray beam. 
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Fig. 4. Determination of the ratio w/L and center channel for a linear detector with 
pixel size of 50 /iin from the variation of the beam position during a scan through 
the primary beam by fitting Eq. 12 to the recorded data. Shown are the channel 
numbers at which the beam is observed vs. the detector angle during the scan. The 
lower plot shows the experimental data as well as the fits when the linear trend is 
subtracted. A fit with (full line) and without (dashed line) considering a detector 
tilt is shown for two detector distances of 250 mm and 380 mm. 
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Fig. 5. Six reciprocal space maps recorded around the Si (331) Bragg peak using 
the same sample movement, however, using different parts of a linear detector for 
detection. The white cross marks the nominal position of the Si (331) Bragg peak. 
Panels (a,d) were recorded when the detector was offset to lower angles, panels (c,f ) 
with an offset to higher angles, whereas panels (b,e) were recorded with the signal 
centered on the detector. For panels (a) to (c) the exact reciprocal space conversion 
described in the text was used, whereas for panels (d) to (f) the "channel per degree" 
approximation was applied leading to errors in the observed peak positions. 
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Fig. 6. Determination of detector parameters. Shown is the optimization error vs. the 
eight considered parameters: the center channels of the detector and the width of 
the pixels in both directions, the detector tilt (azimuth and tilt angle), the detec- 
tor rotation and the offset of the outer goniometer stage. The found optimum is 
marked by a black circle. All other points correspond to fits which did not reach 
the global minimum, those were produced from various different starting parame- 
ters of the fit. Green lines are guides to the eye to visualize the minimum found in 
all the parameters. The points of the fits are colored to enable the identification of 
correlations. 



Synopsis 

Algorithms for the reciprocal space conversion of linear and area detectors implemented in an 
open-source Python package. 
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